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Abstract: Some models for developed turbulence are considered; they are shown to obey a 
large fluctuations theorem, and one among them also obeys a response reciprocity relation of 
Onsager's type. This illustrates and extends ideas and techniques developed in earlier works 
mainly for non equilibrium problems in statistical mechanics. 



§1. Introduction. Earlier results. 

The question: "what are the measures describing turbulence?" has been repeatedly raised 
in a clear form by Ruelle, e.g. see [Rl], p. 6^-8, who proposed on many occasions that the 
probability distributions describing turbulence should share some selected properties among 
the many ones enjoyed by the SRB distributions for axiom A attractors: and in fact it has 
become customary to call SRB distributions the probability distributions on phase space that 
describe the statistics of turbulent (chaotic) motions, [ER]. 

It does not seem that Ruelle ever wrote explicitly how the above prescription could be actually 
implemented for testing: but his idea was made quite clear through his writings and seminars. 
Nevertheless the generality and the breadth of the proposal were never really picked up with 
the purpose of obtaining concrete predictions. 

Recently the following interpretation of Ruelle's prescription for describing turbulence and, 
more generally, motions having an empirically chaotic nature has been proposed, [GC1],[GC2]: 

Chaotic hypothesis: A reversible many particle system in a stationary state can be regarded as 
a transitive Anosov system for the purpose of computing macroscopic properties. 

See [AA],[S],[R2],[Bo] for the notion and properties of Anosov systems. In the very common 
case of non reversible systems one has to replace the transitive Anosov propery with the more 
general Axiom A system property, [R2]. 

The time evolution is described by a flow Vt generated by differential equations on a continuous 
phase space T . 

However we shall always regard the time evolution, or the dynamics, as a map S acting on the 
"phase space" C of the observed events, also called timing events, which could be, for instance, 
the occurrence of a microscopic binary "collision" in a particle system or the event in which a 
prefixed component of the velocity field assumes a certain value, or a maximum value, in fluid 
motions (e.g. a typical example is Lorenz' choice for timing the observations via the maxima of 
the variable he was calling Z, [L]). If x £ C then Vtx evolves in T until the next timing event: 
Sx = V t ( x )X £ C, ift(x) is the time elapsing between the timing event x and the next. 
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The above chaotic hypothesis implies that the macroscopic time averages of observables are 
described by a probability distribution /j, on the "phase space" C of the observed events, also 
called timing events, (which could be, for instance, the occurrence of a microscopic binary 
"collision"). 

The existence of the distribution /j, is assumed a priori in general, as stated by the following 
(extension, see [GC2]) of the zeroth law, [UF], giving a global property of the motions generated 
by initial data chosen randomly with distribution /io proportional to the volume measure on C: 

Extended zero-th law: A dynamical system (C,S) describing a many particle system (or a 
continuum such as a fluid) describes motions that admit a statistics /j, in the sense that, given 
any (smooth) macroscopic observable F defined on the points x of the phase space C, the time 
average of F exists for all jio -randomly -chosen initial data x and is given by: 



where /j, is a S-mvariant probability distribution on C. 

If one assumes the chaotic hypothesis, then it follows that the zeroth law holds, [S,Bo,R2]; 
however it is convenient to regard the two statements as distinct because the hypothesis we make 
is "only" that one can suppose that the system is Anosov for "practical purposes": this leaves 
the possibility that it is not strictly speaking such and some corrections may be needed on the 
predictions obtained by using the hypothesis. The corrections are supposed to become neglegible 
in the thermodynamic limit (in statistical mechanics) or in the large Reynolds numbers regimes 
(in fluids): but the hypothesis, in our intentions, has a meaning similar to that of the ergodic 
hypothesis. The latter too is supposed to be very reliable (under general circumstances as there 
are well known exceptions, e.g. the free gas in statistical mechanics or some very special Navier 
Stokes flows, [M]) in the mentioned limiting situations: but it is very often as reliable in systems 
with few degrees of freedom. 

Still the above principle may look at first sight quite abstract and, in the end, useless. Our 
idea in [GC1],[GC2] was that instead it could eventually reveal itself comparably powerful to 
the ergodic hypothesis in equilibrium non dissipative systems: we began investigating whether 
this statement could be verified in some concrete cases of non equilibrium systems governed by 
a reversible dissipation mechanism and the following results can be obtained. 

(1) The fluctuation theorem (in reversible dissipative systems): The energy dissipation or more 
precisely the entropy creation per timing event <r(x) in the configuration x, identified with 
the logarithm of the inverse of the phase space volume contraction e~ a ^ x \ has a (future) time 
average positive by assumption, but it is a fluctuating variable in the models considered. 

Therefore one can consider the average (a) T (x) of <r(x) evaluated on a time interval of r units 
during which the system evolves between S~ T l 2 x and S T l 2 x: 



Note that the timing events occur at varying time intervals: the time between the event x £ C 
and Sx, denoted above t(x), has an average value to = (t(-)) and the entropy generation rate 
per unit physical (i.e. continuous) time is t$ 1 a(x). 
Reversibility is not contradictory with the presence of dissipation: see [PH], [ECM1], [CG2], 
and the models considered in the present work. 




T-l 



(1.1) 




r/2-l 



(1.2) 



k = -r/2 



2 



We can study the probability distribution of this random variable, with respect to the SRB 
distribution of the motion, i.e. in the stationary state. If we introduce, for convenience, a 
dimensionless fluctuation variable p by setting: 

(7 T (x) = (<r)+p (1.3) 

where p will have (therefore) SRB average 1. 

Then the main result of [ECM2],[GC1],[GC2] is a general property of the probability distribution 

of p denoted ir T (p) = e~ T< * T ( p \ with respect to the stationary statitics \i (i.e. with respect to 
the SRB distribution): 

lim log WT ( p \ = 1 for all p (1.4) 

^m|)t(it)| ir T (—p) 

The above property has no free parameters and it can be tested in some numerical experiments: 
in fact it was discovered in the experiment on particle systems in [ECM2] inspired by the liter- 
ature on the SRB distributions; in [ECM2] the connection with the Ruelle's proposals was also, 
sketchily, worked out. The detailed discussion of the connection between the chaotic hypothesis 
and the fluctuation theorem (1.4) is presented in [GC1],[GC2], with some mathematical details 
worked out in [G1],[G2]. 

The above result applies as well to certain fluid mechanics flows proposed first in [ZJ] and 
mentioned in [GC1],[GC2] and discussed again below, for completeness. The consideration of 
this kind of models goes back to older ideas, see [PH] for a review, which inspired also [ECM2], 
[GC1], [GC2]. 

(2) Onsager reciprocity: suppose that the system is subject to external thermodynamic forces, 
i.e. its equations of motion depend on parameters a, b, c... such that the system is hamiltonian 
when the parameters values are and becomes dissipative but still reversible, when they are 
different from 0. The response current, J , to variations of the parameters is defined together 
with its SRB average j as: 

1 da(x) 

Ja(X) = 7, , Ja = {J a) SRB (1-5) 

to oa 

where to is the average time between timing events. 

Then, if the total energy of the system is kept constant, the reciprocity relation is: 

dbj a \ ab=0 d = L ab = L ba = J d a j b \ ab=Q (1.6) 

and it can be obtained, from the above chaotic hypothesis, in the case of various statistical 
mechanics models, [CG]. For a general discussion of Onsager's reciprocity see [DGM], [ELS]. 
Result (1) is a large deviation theorem and (therefore) it can be tested only in systems with 
few degrees of freedom (because ( T is proportional to the number of degrees of freedom). It is, 
however, a prediction. 

Result (2) is, instead, an a priori test of the hypothesis: in fact reciprocity is independently 
known to hold (for macroscopic systems). It is, therefore, not a prediction but a check if 
consistency with the body of the results that are considered independently established, see 
[DGM]. 

In this paper the above large deviation theorem (already derived as the main result in [GC1], 
[GC2]) is reproduced in the present context, for completeness as well as to introduce the basic 
ideas and notations associated with the chaotic hypothesis. Then we investigate the extension 
of the result (2) to another model related to the incompressible Navier Stokes equation in the 
Kolmogorov Obuchov scaling regime. 
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§2. A model for developed turbulence. 



We consider the Navier Stokes (NS) equations in a box [—I 1 ,! 1 ] 3 , with periodic boundary 
conditions and for an incompressible fluid. If the velocity field is written in a Fourier series as: 

kji 

with 7 complex vectors with 7 = ~y~ (reality of the velocity field) and 7 j. _L k (incom- 
pressibility) then (it is well known and easy to check that) the NS equations become: 

l k = ~ i E (it, • ^2)Hi 7 ki +Rg k -vk 2 l k _ (2.2) 

where n j. is the orthogonal projection over the plane orthogonal to k; v is the kinematic 
viscosity and Rg is the forcing (of course orthogonal to k) which will be taken to be non 

zero only for a few components with small k_ . Since k_ = n with n integer this means that 
the force acts only on the high length scale components. 
For simplicity we may think that the forcing has only two non vanishing components Rg fc0 , 

Rg , corresponding to two linearly independent wave numbers k®, k_° 2 . The simpler case of 

only one non zero component can be trivial {e.g. in 2 dimensional NS when the forcing acts 
on the smallest k, \ k\ = ko, [M]), and is therefore discarded here in favor of the next to the 
simplest. We shall keep g fixed throughout the analysis and we set g = max | g \. 

The number R therefore determines the forcing strength and will be identified with the 
Reynolds number (we also keep the container size L and the viscosity v fixed). We take 
7 = since it is the conserved center of mass velocity. 

In order to obtain equations in the framework of this paper from the phenomenological theory 
of Kolmogorov-Obuchov, [LL], we shall assume that the above equations can be replaced by 
the following simpler ones: 



j k =-i (lk ' i2) n t7t 1 + 9_ k \k\<k f 



a lk~ i 



(2.3) 

• i 2 )n i 7i 1 **<|*|<*ii + CR£</) 1/2 i'" 1 

ti+t 2 =t 1 

where, if ko = ^r, the wave vector kR is the Kolmogorov momentum scale kR = koR 3 ^ 4 , ([LL] 
p. 122, (32.6)). We also set Nr to be the number of wave vectors pairs ("modes") k, — k such 

that k R <\k\< k R + {RLgfl^-\ then N R « (^) 2( ^P^ « R 2 ■ 

Then the total number of pairs k_, ~ k of modes considered is N k, R 9 / 4 among which Nr k, R 2 
are viscous modes, while the others are mertial modes. The number of independent components 
of the field 7 is AN (recall also that, taking into account the reality and incompressibility 

conditions, forces the 7 , 7 to have only two linearly independent complex components). 

We call (2.3) the KO model. 

This means that the equations for the amplitudes 7 corresponding to k's in the mertial 

range, ko < | k \ < kR, are "governed" by the reversible Euler equations. In the viscous range, 
\k\ > kR the dissipation phenomena will be idealized by saying that the equations are simply 
such that only the modes k with kR < | k I < kR + (R Lg) 1 ! 2 !^ -1 have a non zero amplitude 
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and evolve in such a way as to keep the total energy constant. This means that the parameter 
a plays the role of an effective thermostat (or viscosity), which has to be chosen so that the 
total energy is constant, i.e. so that JjXlfc I T fc | 2 = ; 

Sfc L k ' 2_ k 

= ~ i- p 

l^\k\>k R IT ^ I 

The phase space contraction rate is therefore the divergence of the r.h.s. of (2.3), i.e. ANROt(x). 
Hence the entropy generation per per timing event is a(x) such that: 

1 E k f h ' T _ j~ f] P f six) 

r w = <■> = (^^,:,h t p - (2 - 5) 

where to is the average time between timing events and e(x) and kT(x)/2 are simply defined 
respectively by the numerator and denominator of the above fraction defining a(x). 
In (2.5) and in the following we neglect, for simplicity, the variability of the time t(x) elapsing 
between timing events. If one wished to avoid this approximation the r.h.s. of (2.5) should be 
replaced by t$ 1 times its integral over the continuous time trajectory over [0, t(x)] of AN Ra(Vtx) 
(using the notations fixed in §1). 

The Kolmogorov length k^ 1 is introduced here phenomenologically and we do not attempt at a 
fundamental derivation of (2.3), (2.4). Therefore (2.3) has to be regarded as a phenomenological 
equation. A class of similar models was introduced in [ZJ]. 

Note that a is proportional to the work s(x), per unit time and per viscous degree of freedom, 
performed on the system and dissipated into heat in order to keep the total energy constant: 
the proportionality constant is 2/kT(x) with ^kT(x) = ^ k \ <kR \ j k | 2 (which, however, 
is not a constant on the motions described by (2. 3), (2. 4) because of the imposed constraint that 
Si IT J 2 is constant, rather than Y,\k\>k R llj 2 )- 

Hence can be thought of as the average amount of energy dissipation per unit time by the 
flow divided by the kinetic energy per mode contained in the viscous modes. The first quantity 
plays a major role in Kolmogorov's theory, see [LL] p. 119, and its average is usually called e, 
(see [LL], (31.1)). Since the kinetic energy contained in the viscous modes can be considered 
as a kind of "temperature" we see that (&) + , is proportional to an entropy "production rate". 
More appropriately we can say that, for R large, (c) + is proportional to the "energy dissipation 
rate" over a kinetic quantity equal to the average kinetic energy contained in the viscous modes 
provided, for large R, the two quantities can be regarded as independent random variables (so 
that the average of this ratio equals the ratio of their averages). 

Note that for the above model (2.3) the time reversal map i : { j } — ► {—7 k } has the 

property that iV-t = V_t if Vt is the flow generated by the (2.3): see §1 for the notation. 

As discussed in §1, we shall call C the space of the velocity fields 7 with a fixed value U for 
the total kinetic energy and such that a fixed velocity component {e.g. an arbitarily selected 
viscous component 7 j. 1 in the ^-direction) has a given value (or a local maximum, to follow 
[L]). We call the space C the space of the timing events as we shall imagine to record the velocity 
field every time one event x in C occurs. Then the time evolution flow induces a map S on C 
mapping one element x £ C into the event of the same type which occurs after it. If t(x) is the 
time between the timing event x and the next following it it is Sx = V t ( x )X. 

The time reversal symmetry for the continuum flow becomes the symmetry Si = iS -1 on the 
timing map S that we shall call the evolution map. 

A similar model is the one considered originally by [ZJ]: in this model the equations of motion 
are given by the (2.3) but with the thermostatting forces —ay present also in the first 
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(2.4) 



equation in (2.3). The value of a is still determined by imposing the constancy of the total 
kinetic energy. Therefore the (2.4) is modified by replacing the sum in the denominator by 
an unrestricted sum. This means that in this model T(x) is rigorously constant and therefore 
there is a rigorous proportionality between the phase space contraction rate and the energy 
dissipation rate. 

§3. Chaotic hypothesis and vanishing Lyapunov exponents. 

In the above models the dimension of the phase space C on which the time evolution is repre- 
sented by a map between timing events is two units less than the dimension of the space of the 
velocity fields (where the motion is described by differential equations like (2.3)). 

This reduction is useful not only because it eliminates a degree of freedom which is "trivial" 
(by fixing the energy which is, in any event, constant) but also because it eliminates the degree 
of freedom corresponding to the direction of the flow which is responsible for the existence of a 
Lyapunov exponent. 

Obviously had we not assumed the point of view of regarding the time evolution as a map on 
C, then this would have caused us the problem that the chaotic hypothesis would have been 
in conflict with a basic property of Anosov systems: in fact the Lyapunov exponents for such 
systems are separated from by a gap A > 0. 

But one may well have doubts that considering the evolution on a space of timing events is 
sufficient to avoid that the chaotic hypothesis comes into a manifest contradiction because of 
the possible existence of other vanishing Lyapunov exponents. 

For instance in some examples, e.g. model 2 in [GC2], the dimension is reduced by 4 units 
because the equations of motion preserve a component of the center of mass momentum, pro- 
vided it has value, and also the corresponding center of mass position coordinate. If the total 
momentum is not initially it relaxes to the value with a negative Lyapunov exponent, while 
the center of mass motion relaxes to a rectilinear motion with zero exponent. 

Therefore fixing the values of the total momentum component and the center of mass coordi- 
nate further reduces the dimension by 2 and eliminates two more "trivial coordinates". 

In this case the conservation of the horizontal momentum is due to the special boundary 
conditions used: it was pointed out in [GC2] that changing the boundary conditions may turn 
the horizontal momentum and the center of mass position into a non exactly conserved quantity, 
no matter which is the value given to it initially. And in such circumstances we can expect that 
while the momentum evolves towards its stationary value (0) at exponential rate (i.e. "with a 
non zero Lyapunov exponent") the slow almost linear motion of the center of mass generates a 
vanishing Lyapunov exponent, at least in the limit of large systems). 

But one should expect that when the system becomes large (i.e. in the thermodynamic limit) 
the boundary conditions should have neglegible influence on the macroscopic properties of the 
system. Therefore the macroscopic dynamics may be insensitive to the above variability and, in 
consequence, adding as a constraint that the horizontal momentum vanishes exactly, hence the 
center of mass position is also on rectilinear motion exactly, should not affect the macroscopic 
behaviour. At least if the forces needed to impose the constraint are provided by a minimal 
constraint principle (see below). 

One should note that the above replacement of an almost conserved quantity by an exactly 
conserved one has in fact the consequence of eliminating a trivial pair of Lyapunov exponents: 
one is and one is negative and describes the relaxation to equilibrium of a macroscopically in- 
teresting quantity (the total horizontal momentum, i.e. the horizontal current, in the mentioned 
example). 

In [GC2] it was proposed, see §8, that the above mechanism might be quite general. In 
the sense that a non equilibrium system may have many, perhaps very many Lyapunov 
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exponents (evidence in this direction can be found in the basic paper [LPR]): this does therefore 
(apparently) violate the chaoticity hypothesis of §1 strongly. 

In such systems the chaotic hypothesis can still be assumed, it was suggested in [GC2] §8, if the 
dynamical variables responsible for the existence of the Lyapunov exponents can be identified 
togheter with the conjugate variables (which would have negative Lyapunov exponents, and 
which are expected to exist quite generally at least if one accepts as generally valid the pairing 
rule discovered in [ECM1] and discussed in [EM], [ES], [GC2]), and then one modifies the 
equations of motion so that the identified variables are exactly constant (with the one with 
negative Lyapunov exponent fixed at its equilibrium value). 

The proposal in [GC2] was that the minimal forces necessary to impose the constraint would 
not affect the macroscopic behaviour of the motions: here by minimal forces we proposed to 
intend the forces prescribed by Gauss principle of least constraint, (see [LA], vol. II2, p. 470, 
and appendix Al below). The idea being that the dynamics enforces the constraint whether 
it is present or not, provided the constrained variables are assigned the appropriate stationary 
average value. 

The chaotic hypothesis can then be used to describe the evolution of the remaining coordinates, 
if one has taken into account all the macroscopic constraints so that no Lyapunov exponents 
are present any more. 

The question is whether one can identify concretely all the dynamical variables that can be 
fixed to have a well defined value without affecting the macroscopic behaviour of the system. 

A guide to finding such variables should be the macroscopic equations that the system is 
supposed to obey. To give an example of what I have in mind consider a dense gas which 
can be regarded as an incompressible inviscid fluid. In this case we can reasonably expect, 
[EMY,KV,LOY], that the particles motion proceeds in such a way that the macroscopic average 
velocity field u(x) evolves according to the Euler equations, while the displacement field, 
relative to the initial configuration, changes accordingly. 

Therefore given the solution of the Euler equations corresponding to the initial macroscopic 
state of the gas: t — ► u(x,t), one can impose that the average velocity locally coincides with 
u(x,t). This means that we can consider a lattice of cubes A 7 with microscopically large and 
macroscopically small side, see [KV], and call x_ the center of the cubes; then we can impose 
the constraints: 

E^(£p|f = «U 7 ,*)> f E^(iP(|f- u) 2 = N,kT (3.1) 
3 3 

where T is the temperature and x-y is the characteristic function of A 7 , via the imposition of 
auxiliary forces derived from Gauss minimal constraint principle (see [LA]), i.e. by writing the 
equations: 

1 3 = ^£,-> Ej = fj + £ U 7 + - « ))x 7 (?p (3-2) 

7 

where a 7 are determined by imposing that the constraints are exactly verified. This is: 

a_ 1 = u{x ,t) (3.3) 

and /3-y is determined likewise. 

One expects that in this way a very large number of coordinates which evolve with a negative 
Lyapunov exponent is eliminated; and that the corresponding position coordinates, i.e. the 
coordinates of the displacements with respect to the initial positions, change with a Lyapunov 
exponent. 
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Some care has to be exercized here: the remaining coordinates will be supposed to vary with 
Lyapunov exponents separated by a gap A from the value 0, independently of the system size. 
This means that the evolution of the very large number of remaining coordinates takes place 
on a very fast time scale A -1 . We can expect the latter to be the same time scale on which 
the local averages reach the "equilibrium" value given by u(x_,t) (this is a consequence of the 
pairing rule, [ECM1],[GC2], if accepted). On the other hand the displacement variables should 
evolve with a Lyapunov exponent: but the latter is only if compared to X. It may well 
be non zero and very small so that its inverse is a macroscopic time, as it is shown in many 
numerical experiments (that show positive Lyapunov exponents of the "lagrangian motions"). 

This would mean that the chaotic hypothesis is valid for finite systems, in the space of the 
timing events, and that it becomes invalid strictly speaking only in the thermodynamic limit, 
when however it is not valid for trivial reasons and assuming it simply corresponds to think that 
some approximate conservation laws have become exact macroscopic evolution laws thus imply- 
ing the vanishing of some Lyapunov exponents. The phenomenon can be avoided by imposing, 
already in finite systems, the approximate conservation laws as exact laws and eliminating the 
relative coordinates (thus restoring a uniform gap in the Lyapunov spectrum). 

In our fluid model the situation is likely to be similar: and one has to interpret accordingly 
the chaotic hypothesis: namely one has to think that there are no Lyapunov exponents or 
that if they are present they can be eliminated by adding extra forces that fix the value of 
some observables which would relax slowly to equilibrium (thus generating vanishing Lyapunov 
exponents) without affecting the behaviour of the system (except of course for what concerns 
the long time correlations of the very same observables evolving with vanishing Lypapunov 
exponents). 

As an example of the above discussion one can argue that the model (2.3) and its modification 
in which the thermostatting forces act on all modes (considered at the end of §2, [ZJ]) can 
be considered without obvious contradictions to verify the chaotic hypothesis for finite R, 
and uniformly in R at least as far as the properties of the observables that relax quickly to 
equilibrium are concerned. 

§4. The SRB distribution. 

The application of the chaotic hypothsis, as we proposed in [GC1],[GC2], is very similar to 
the applications of the ensembles method in equilibrium statistical mechanics. In that case 
one does not really need to evaluate all the monstrous integrals over phase space to deduce 
remarkable macroscopic properties: a great example is provided by Boltzmann's derivation of 
the heat theorem, [B], [G3], i.e. of macroscopic thermodynamics. 

The chaotic hypothesis can perhaps be used in a similar way because it provides us (see 
below) with an apparently impossibly complicated expression for the SRB distribution: it is 
nevertheless an expression on which, as we have shown in [GC1],[GC2], one can work quite in 
detail and from which the concretely testable propeties discussed in §1 have been derived. 

The key point is that the hypothesis implies the possibility of defining a natural coarse graining 
of the phase space which is also mathematically precise (a problem that has been previously 
repeatedly debated without ever reaching clear conclusions, in my understanding at least, [G4]). 

In fact in an transitive Anosov (or axiom A) system one can define a partition £ of phase space 
C into Af sets E\, E2, . . . , Ej\f which are "regularly" shaped (i.e. with non empty interior, with 
boundary which has zero volume and which is smooth at least in the sense of Holder continity). 

The sets in £ have the (highly non trivial) property that if one defines a compatibility matrix 
C'ij by setting C'ij = 1 if the interior of SEi intersects the interior of Ej and otherwise, then 
any sequence j = (jk)k=-oo,oo such that Cj k j k+1 = 1 for all k's (compatible sequence) is the 
history on £ of one and only one point x £ C, in the sense that S k x £ Ej k for all k's. 
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Viceversa each point x £ C is generated by one compatible sequence, naturally called the 
history of x on £, in the sense that it is the only point in the intersection C\kS~ k Ej k . There 
may be more sequences determining the same point but this happens for a zero volume set of 
points (this is a quite trivial ambiguity as it is similar (and in fact very closely related to) the 
well known ambiguity that one has in representing the reals in basis 10 by digits that end in 
a sequence of 9's or of 0's). Furthermore the matrix C admits a power C q with all the entries 
positive (this is a consequence of the transitivity). 

The above property allows us to think of the compatible sequences as a representation of our 
phase space and to regard the volume measure as a measure on the space of the compatible 
sequences. The SRB distribution \i also can be regarded as a probability distribution on the 
space of the compatible sequences. The latter, in turn, can be conveniently regarded as the 
space of the spin configurations of a one dimensional spin model (the spin at the site k being 
jk)- 

The distribution \i has a remarkable representation when considered as a distribution on the 
space of the compatible sequences, 

To describe the representation we need some further notations. Given a symbol j we can find a 
semiinfinite compatible sequence j = (j, j' , j" , . . .) whose entries depend only on the value of j; 
likewise we can find a seminfinite compatible sequence j =(..., j" , j' , j) whose entries depend 
only on j. This poperty, consequence of the above mentioned transitivity, can be used to locate 
conveniently a point whose symbolic representation contains a string j a]j = ja, ja+i, ■ ■ ■ , jb- one 
simply continues the string j ai into the infinite string: j = ( j (j a ), j a b , j + (jb)) obtained 
by merging sequentially the three strings j (j a ), j ab and j + (jb)- 

Let E = Ej_ Mt ^j M = Pi%L_ M S~ k Ej k : this is the set of points x such that S k x £ Ej k for 
k = j-M, ■ ■ ■ ,3m- We fix a point xe = %i- M ,...,j M £ Ej_ My ...j M to be the point whose history 
is obtained by continuing arbitrarily the sequence j-M, ■ ■ - ,3m "to the right and to the left" 
into a biinfinite compatible sequence j as described in the previous paragraph. The point xe 
will be called the center of E. 

The expansivity of an Anosov map implies that the sets E = Ej_ My ^j M are very small (their 
diameter tends to as e~ XM if A is, essentially, the gap in the Lyapunov spectrum) and we can 
define a distribution hm by assigning a weight to each of such sets. 

The weight that we assign to them is related to the expansion coefficient A Ui m(%e) of the map 
S M as a map between S~ m I 2 xe and S m I 2 xe- 

The expansion (contraction) coefficient A u (x) (respectively A s (x)) of S, at x, is the jacobian 
determinant (evaluated at x) of the transformation S as a map of the unstable (stable) manifold 
W u (x) (W s (x)) into itself. For a discussion of the notion of stable and unstable manifolds see 
[Rl], [ER]. Therefore the expansion (contraction) coefficient A Ui m(%) (respectively A s m(*)) 
are given by: 

M/2-1 _ M/2-1 _ 

VW= II MS^) = e Mj ^ x \ A S) mW= A,(S>z) =V Ml -« (4.1) 

j = -M/2 j = -M/2 

and the distribution hm is defined by giving to each Ej £ £m, with center xe 3 =%j, a 
weight proportional to the product of the expansion coefficient A u m( x ]) times the inverse 
of the sine of the angle {)(S M l 2 Xj) formed by the stable and unstable manifolds at S M l 2 Xj\ 
^m(*j) d — am , d(S M l 2 Xj). So that the integral of a smooth function F is: 

(a \t?( \ de f A ",MV 1 ( a; i) F ( x i) ( , 

H M {ax)F{x) = — 1 — (4.2) 

C J2j ^u,M( x j)^M ( x j) 
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This is interesting because the following theorem holds: 

Theorem: If (C, S) is a transitive Anosov system the SRB distribution /j, exists and the /j, average 
of a local function F is: 

/ n{dx)F{x) = lim / /i M (dx)F(x) (4.3) 

Jc M^co J c 

and "local" means that F depends "exponentially little" on the symbols with large time label in 
the symbolic representation of the phase space points (see above; i.e. \F(x) — F(y)\ is exponen- 
tially small as k — ► oo if the histories of x,y at the times between —k and k). 

The above is a trivial reformulation of a deep result of Sinai. It was pointed out in [G1],[G2] 
and used in [GC1], [GC2]. I will call it Sinai's theorem. 

The original statement is that the SRB distribution /j, exists and it is a Gibbs state with 
potential log A" 1 (a;): see [Bo],[R2],[S] for a discussion of this form of the statement. 

The connection between [S] and the above formulation is discussed in [Gl] where (4.3) is 
discussed with \im defined as in (4.2) without the factors 8 M 1 (xj). In spite of the apparently 
strong modification the extra factor 8 M 1 (xj). introduces, it is easily seen that (4.3) is valid by 
examining the proof of (4.3), see [Gl], [G2]. 

The proof is based on the interpretation of (4.2) as a probability distribution on the space of the 
compatible strings. In this interpretation one immediately recognizes that (4.2) corresponds to 
a finite volume Gibbs distribution for a suitable short range hamiltonian defined on the space 
of compatible strings. The extra factor 8m (xj) corresponds to considering the same Gibbs 
distribution just with a different boundary condition: this becomes irrelevant in the limit as 
M —> oo because one dimensional Gibbs states with short range interactions do not have phase 
transitions and therefore are i insensitive to changes in the boundary conditions. 

The choice (4.2) as an approximating distribution for /j, is much better than the one without 
the factors 8m{xj)~ 1 because it leads to simpler formulae and arguments: we shall call (4.2) 
a balanced approximation to the SRB distribution because as we shall see it is reminiscent of 
a probability distribution verifying the detailed balance (which however is not verified in our 
models, except in forcing, i.e. in equilibrium). 

In the case in which the system is reversible (as the model in §2), i.e. when there is an isometric 
map i : C — ► C such that iS = S~ 1 i and i 2 = 1, one can add to all the above properties the 
further properties that £ can be taken time reversible. This means that if E £ £m then also 
iE £ £m, with i being the time reversal operation, and ixe = XiE- Furthermore, recalling 
(1.2), the following symmetry holds: 

a M {x) = -a M (ix), A UtM (ix) = A s ^ M (x), 8 (ix) = 8 (x), 8 M {ix) = 8_ M {x) (4.4) 

which are identities (see [GC2]) simply implied by the definition of a or by the fact that the i 
operation changes the stable manifold for S at x into the unstable manifold at ix (still for for 

Furthermore the following relation holds between A u , A s , 8, a: 

A u (x)A s (x)^ = e~«*\ A UiM (x)A SiM (*)P^\ = e~ M ^ = e~ M ^ (4.5) 
8{x) 8-m\x) 

if one recalls the definitions (1.2), (1.3); the second relation follows from the first: the l.h.s. is 
in fact the phase space contraction under the map S M evaluated at the point S~ M l 2 x. 
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For small dimensionality systems the sets E q can be quite easily and effectively constructed: 
but the construction, which involves among other things, solving the equations of motion for 
many data, becomes quickly practically impossible, [FZ]. 

§5. Applications of the Chaotic hypothesis. Fluctuation theorem for the KO 
model. 

We study the KO model in (2.3), or the modification considered at the end of §2, and we 
suppose that the system is kept in a stationary state at a constant energy U under the action 
of a force field g and of the thermostatting mechanism provided by the terms —ay in the 
equations of motion. 

We study the probability that the fluctuaton variable p is in a small interval I q = [q,q + dq\. 
We use the notations and the approximation \i T to \i described at the end of §4 (see (4.2)) with 
F(x) = <T T (x); the probability that p £ I q over the probability that p £ I_ q is, for large r: 

Since \i T in (4.2) is only an approximation at fixed r an error is involved in using (5.1). It can 
be shown that this error can be estimated to affect the result only by a factor bounded above 
and below uniformly in T,p, [GC1], [Gl], [G2]. This is a remark technically based on the proof 
of the theorem quoted in §4 (which relates the problem to the properties of a one dimensional 
short range Ising chain, a technical tool that is usually called the "thermodynamic formalism") 
and it is valid in general for any system verifying the chaotic hypothesis, i.e. for any reversible 
transitive Anosov diffeomorphism, [G2]. 

It is now possible by using the reversibility to establish a one to one correspondence between 
the addends in the numerator of (5.1) and the ones in the denominator, (aiming at showing 
that corresponding addends have a constant ratio which will, therefore, be the value of the ratio 
in (5.1)). 

The ratio (5.1) can be written simply as: 

ZE^pei-JuU^rHxj) ~ E B , pe / 5 A.,,fe)«:J(^) 

where xj £ Ej is a point in Ej . In deducing the second relation we make us of the existence of 
the time reversal symmetry i and of (4.4). 

It follows then that the ratios between corresponding terms in the ratio (5.2) is equal to 

A u \.(x)A s ^(a;) ^-l^ 1 '^ ■ This is, by (4.5), the reciprocal of the total variation of phase space 

volume over the r time steps, if the evolution is regarded as a map on C, between the point 
S~ T l 2 x and S T l 2 x: but the reciprocal of the total phase space volume contraction over a time 
r is e q ^ + T . Hence the ratio (5.1) will be e T ^+ q proving (1.4). 

It is important to note that there is one error ignored here. As pointed out in the discussion 
above the use of \i T to evaluate the probability is not immediately justified by the theorem of 
§4 as the function F of which we study the distribution is not a "local" function, because cr T (x) 
in (1.2) is not localized. In fact it depends on the history of x for a time between — r, r (hence 
we essentially need the knowledge of the symbols between — r/2 and r/2 in the history of x 
to compute the value of a T (x)). This is a delicate point which has been discussed in detail in 
[G1],[G2]: showing that this interchange of limits problem is really not a problem at all requires 
going into the details of the proof of the theorem of §4. 
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(5.1) 



(5.2) 



One should note that other errors may arise because of the approximate validity of our main 
chaotic assumption (which states that things go "as if the system was Anosov): they may 
depend on R, i.e. on the number of degrees of freedom, and we do not control them except for 
the fact that, if present, their relative value should tend to as R — ► oo. 

The p independence of the limit in (1.4) is therefore a key test of the theory (for a Anoson 
system; and from the detailed estimates in [G2] one sees that the limit is reached as r — ► oo 
with corrections of order 0(r _1 ), for p in a fixed bounded interval). 

§6. Reciprocity in a model for a shell motion. 

In [CG] it will be shown that the chaotic hypothesis implies also, in a variety of non equilibrium 
statistical mechanics reversible models, the Onsager reciprocity relations. 

It is tempting to apply the above ideas to the models, so far considered, for fluids in states 
of developed turbulence. There is, however, a basic difficulty: namely the above mentioned 
derivations deal with "infinitesimal deviations from equilibrium": because they express proper- 
ties of the second derivatives of the dissipation rate with respect to the thermodynamic forces 
evaluated at thermodynamic forces. 

The fluid models considered above (KO and its variation mentioned in §2) do not fulfill the 
condition of being close to equilibrium: in fact we always imagine R to be very large so that 
the Kolmogorov Obuchov theory applies and the models may be regarded as describing chaotic 
motions and as physically reasonable. 

There is no known extension of Onsager reciprocity to strongly forced systems, [DGM]: there- 
fore we shall not insist in studying the above models. But a possible application to fluid motions 
can be found by considering a related model. 

We consider a fluid whose velocity field contains only components with momenta in the range 
M„ defined by 2"~ 1 ko < | k | < 2"ko containing N pairs k, — k: this is often called a momen- 
tum shell. The equations of motion will be: 

11=-*' E (l il ^2)n i 7t 1 + / i -a7 i 2 n ~ 1 k < \k\,\k 1 \,\k 2 \)<2 n k (6.1) 
ti + t 2 =t 1 

and a is fixed so that the total kinetic energy is 2NkT, hence the relations between kinetic 
energy, the "friction" a, and the phase space contraction exponent a are: 

-Y\j , I 2 = 2NkT, a=— = a= ~ , -, ^t (6.2) 

2^'-i' ' 2NkT ' kT/2 v ' 

if to is the average time interval between timing events. As already mentioned after (2.5), 
for simplicity of exposition, we neglect that there is a difference between the actual time t(x) 
elapsing between the timing event i£C and the successive Sx and that during such intervals 
cr(Vtx) is not strictly constant. In fact a should be really written as an integral over the 
continuous time trajectory t — ► Vtx described in the time t(x) while the point x evolves between 
x and Sx. 

One can think that this is a crude model for the motion of the velocity field components with 
momenta of scale n, if the range M n is in the inertial range of a strongly turbulent flow. 

In this model we can ask the question of the reciprocity of the fluctuation responses with 
respect to the variations of the forces / . 

Denoting /3 a pair k_ , h with h = 1,2,3, and setting fp = fk,h and dp = d ^ , we want to 
show, under an extra assumption (see the two paragraphs following (6.14) below), that the 
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chaoticity hypothesis implies: 



dp> (dpa)^ 



def 



f = o 



Lp-,p d = dp (d'pa)^ 



f = o 



(6.3) 



where /j, denotes the SRB distribution. 

In spite of its appearance (6.3) requires some discussion on its meaning. 

The function t^ 1 a(x) is, by definition, the phase space entropy production rate, it is naturally 
defined on the phase space C. However if > (as it is the case for / ^ ) the motion 

is dissipative and the attractor for the motion has zero volume, see [GC2]. In fact we expect, 
believing the pairing rule discovered in [ECM1], see also [GC2], that the fractal dimension of 
A is macroscopically lower than that of C. 

Therefore the "relevant" values of the function a(x) are those for which x £ A. But A as well 
as the SRB distribution /j, are /-dependent. Hence if we want to discuss the / dependence of 
a we must think that a is defined on a surface that generates the attractor, e.g. the unstable 
manifold of a fixed point O (or periodic motion) W u (0), and which depends on f . 

The consequence is that dp(a(x)) cannot be identified with the partial derivative of (6.2) 
(i.e. jTpJp) but one has also a contribution d "g^ ■ -§j^ as the phase space point x has to be 
fixed on the attractor and it changes therefore with / . 

A simple way to see this from a different point of view is to think of the attractor as represented 
via the symbolic strings j associated with a Markov partition: x = x(j), as discussed in §4. 

The SRB distribution becomes a probability distribution over the family of compatible strings 
j and, in fact, it is a Gibbs state corresponding to a potential which has short range (i.e. by 

a potential with energy per site given by the function logA^ 1 (x(j ))). 

In this language the dynamics becomes f -independent and a becomes a function of the string 
a = a(x(j)). Therefore a depends on / for two reasons: 

(1) a depends "explicitly" on /, see (6.2). 

(2) x( j) depends on / as well, because the correspondence between strings j and points x( j) 
depends on the dynamics. 

Hence we shall use the more appropriate notation: a(x) = s( f ,x( j )) and: 



da ds ds dx(j) 



Of course it is s(0,x(j)) = so that we get: 



da 



f = o 



ds 



f = o 



(6.4) 



(6.5) 



The proper interpretation of (6.3) is obtained by defining dp as in (6.4), which makes a difference 
even when we may only be interested in evaluating Lppi at / = 0. In fact (6.3) involves the 
second derivative of a, for which a relation like (6.5) does not necessarily hold. 
We shall see below that (6.3) holds with the latter interpretation of dpa . Hence we argue that 
the correct definition of the current jp when / 7^ seems to be: 



1 



3P 



(6.6) 



where dpa is defined as in (6.4). 
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One could also object that since the code j — ► x(j) depends on / it might happen that 
the space of the compatible sequences itself changes with /: i.e. that the compatibility matrix 
depends on / introducing a further dependence of a on x and in fact a dependence on / on the 
symbolic dynamics itself. Although this might indeed happen at "large" (perhaps moderately 
large) values of / it is a consequence of Anosov structural stability theorem, see [AA], that the 
compatibility matrix C for the compatible strings j does not change for small enough variations 
of the dynamics, hence of / . In fact establishing the constancy of the stability matrix is the 
first step in the proof of the structural stability (and essentially the reason for the stability 
itself). 

The above definition (6.6) of the notion of current associated with a thermodynamic force 
seems to be new and it should be carefully understood and tested in the nonlinear regime. Here 
we assume by definition that the current associated with the thermodynamic force component 
fp by (6.6) with the derivative of a defined by (6.4). 

The proof of (6.3) now follows a natural scheme, see also [CG]. We set, in an effort to simplify 
notations in the coming formulae, l u>T , l s>T as: 

I u ,t(x) = log A u 1 T (x)6~ 1 (x), l StT (x) = log A StT (x)6zl(x) (6.7) 

so that by (4.4): 

I u ,t(x) - l SlT {x) = Tcr T (x), l u , T {ix) = l SlT {x) + Ta T (x) (6.8) 

and we see that jp = -^(dpa) (see (6.4)) is the limit of the r.h.s. of (4.2) with F(x) = dpa(x) 
and M = r as r — ► oo. We can also say that (dpa) is the limit of (4.2) with F(x) = dp<r T {x), 
replacing a by its average, see (1.2), between —\t and -|t. This is easily justified if one recalls 
the symbolic dynamics interpretation of \i T (discussed after the theorem in §4) in terms of a 
1-dimensional Gibbs distribution for a short range Ising model (see the second comment after 
(4.3)). I will not discuss this (minor) technical point further. Hence: 

Ej A tl ,r(*i)^ 1 ( ;C i)( 5 /3'/3' T r(*i) + p> l u tT {x j)8 p(T T {x j) 

lim todpijp = =^ ■ 

Ej K,r( X j) 6 r (Xj) 
_ Ej K.l^jW^^^dp'lurjxj) Ej K.l^jW^^^dp^rjxj) 
Ej K.kxj^r^Xj) J^j A u, 1 r( X i) 6 r 1 (x j ) 

= lim Udp'p<T T ) + ((dp<l UiT (xj)dp<T T ) - (dp<l UiT )(dp<T T )) 



(6.9) 



By using the time reversal invariance we see that: 

7^ A,, l6~ 1 (xi)dB'L l T dea T 
(dp-l UiT dpa T ) = ^ U ' T T 1 3 J {i U ' T {i = (6.10) 

_ Ej (K,l( x iK\ x j) d l3' l u,r(x j )dpa T (x/) + A u l(ix j )6-\ix j )dpd UtT ^^^ 

where Z denotes the "partition sum" , i.e. the sum in the denominator of (6.9), and the averages 
are with respect to the distribution h t /2,t- Recalling (6.7) this becomes: 

Ej (K.l^iK^^jWli'^.rixj) - As^x^Szlix^dpds^xj)) dpa T (xj) 

~z (6 - n) 
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Equation (6.8) permits us to reconstruct rdpia T ou of the two addends in (6.11). And the 
derivatives at / = can be computed immediately by noting that in such case, by (4.5), 
(6.7),(6.8): 



{dp>l u ,T df}cr T ) - (dp>l UtT )(dpa T ) 



f = o 2 



(dp>a T df}cr T ) - (dp>a T )(dpa T ) 



f = o 



(6.12) 



The quantity (dppicr T ) will converge to a limit (di,picr) because the SRB distribution is stationary, 
neglecting an exchange of limit (see comments below); then by using (1.2): 



dpi ip = lim {dppi a) 

F 1=0 r^oo \t 
t/2-1 t/2-1 

— EE {(df^{S m -)dp<r{S n -))-(dp,<r{S m -))(dp<r{S n -))) 



m = — r/2 n = — r/2 

1 1 00 

T^' a K n+^T E ((dp-<T(S m -)dp*(-))-(dp,*(-))(dp*(-))) 



m = — oo 



f = o 



f = o 



(6.13) 



where the averages in the r.h.s. are with respect to /j, = lim//,-. This, again apart from a 
problem of interchange of limits (see comments below) becomes: 



dp'jp 



1 CO 

f =W E ((^<7(5 m .)^<T(-))-(^<T(-)>(^<T(-)>: 



1 1 °° 

r^" 7 > = ^r E <^^ m -)9^(-)> 

tn Zip, 



(6.14) 



m = — oo 



where the averages are with respect to the SRB distribution (i.e. with respect to the limit of 
fi T ) and the missing terms vanish because of time reversal symmetry (for instance j^(dpp'u) is 
seen, considering (6.2) and (6.4), to coincide with the expectations (dpjpi +dpijp), linear in 
7 and therefore vanishing by the time reversal symmetry). 

The problems of interchange of limits are easily solved: under our assumption that the system is 
a transitive Anosov system the correlations of smooth observables decay exponentially (because 
they become local observables in the symbolic dynamics interpretation of the evolution, provided 
by the Markov partitions), not only for /j, but also for \i T (in the natural sense in which this 
may be interpreted in a finite r case; e.g. by regarding the interval [— -|, -|] as a circle), and 
uniformly in r. 

Note that here there is one hidden assumption: namely the "local observable" to which we want 
to apply the above argument is dp<r which, by (6.8) contains besides the really local part given 
by the first term in the r.h.s. a second part containing dpx( j ) (and dpx(-d k j ) with «9 being the 
shift map on the symbolic sequences): therefore we must assume that also dpx( j ) is "local" in 
the sense of the theorem in §4, which seems a reasonable assumption. One can conjecture that 
it holds in general for any transitive Anosov system (or axiom A system) depending smoothly 
on a parameter fp . 

This also shows that (6.14) is indeed correct and, therefore, recalling that Lppi = tndpijp we 
get (6.3) as well as the important property that the matrix Lppi is positive definite (see the last 
of (6.14)). And one can check (by (6.2) and (6.4)), that the relation Lppi = Lpip is, in fact: 



dp-(jp) 



f = o 



dp(jp-) 



f = o 



(6.15) 
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which is easier to interpret. 
Remarks: 

(1) Note that (6.13) seems to have one "extra term", ^(8/3/3* &) , with respect to what one 
usually expects from linear response theory: this is due to the fractality of the attractor and it 
arises because od the "extra term" in (6.4). 

(2) If one thinks that the above is a reasonable model for the evolution of the velocity com- 
ponents realtive to a shell in the inertial range, then above relation should be subjected to an 
experimental test. In any event the relations are very likely testable in numerical experiments 
on the solutions of the equations (6.1). 

(3) One can remark that if jp =jk,h f° r ft = (k,h) the above reciprocity relations can be 
written as the symmetry of the matrix dp (7/3')/"| v -o on ^ ^ we consl der transformations of 
the system which vary / but keep the total energy U = NkT constant. 

(4) One certainly wants to consider also transformations in which U changes according to 
some equation of state U = </(/)• The function g depends on the physical situations, i.e. on 
the actual mechanism of dissipation that generates the model (6.3): however if the model is 
regarded only as a mathematical fiction, then of course the equation of state is arbitrary. In 
this case, given g, the above derivation of the reciprocity relations still holds up to some natural 
modifications, but the symmetry 8/31(8/3(1) does not directly mean the symmetry (6.15). One 
can, however, check that it implies at / = : 

dfAilP - = 8 f 3{( lf 3, - ^d p U)) (6.16) 

which still implies (6.15) (again by the time reversal symmetry). This is in agreement with the 
philosophy exposed in §3. 
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